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ABSTRACT 

To disdosc detailed genetic mechanisms in hepatocellular carcinoma 
CHCO with a view toward development of novel therapeutic targets, we 
analyzed expression profiles of 20 primary HCCs and their corresponding 
noncancerous tissues by means of cDNA microarrays consisting of 23,040 
genes. Up-regulation of mitosiS'promoting genes was observed in the 
minority of the tumon examined. Some genes showed expression patterns 
in hepatitis B virus-positive HCCs that were different from those in 
hepatitis C virus-positive HCCs; most of them encoded enzymes thai 
metabolize carcinogens and/or anticancer agents. Furthermore, we iden- 
tified a number of genes associated with malignant histological type or 
invasive phenotype. Accumulation of such data will make it possible to 
define the nature of individual tumors, to provide clues for identifying new 
therapeutic targets, and ultimately to optimize treatment of each patient. 



INTRODUCTION 

Primary HCC"* is one of the most common malignancies in the 
wortd. Despite development of novel therapeutic methods in recent 
years, prognosis of advanced HCC remains very poor. Major risk 
factors for HCC are chronic hepatitis resulting from infection with 
HBV or KCV, and exposure to various exogenous carcinogens in- 
cluding anatoxin Bl (1). Molecular approaches have recently revealed 
involvement of altered TP5S, CTNNBI (P-catenin). and/or AXINI 
genes in hepatocarcinogenesis (2. 3). However, these genetic changes 
do rK)l precisely reflect the biological nature of cancer cells or the 
clinical characteristics of individual HCC patients. Like other cancers. 
HCCs marufest diverse clinicopathological and biological phenotypes 
including grade of differentiation, proliferation rate, ability to invade 
vessels, potential for metastasis, sensitivity to chemotherapeutic 
agents, and so on. Hence, analysis of expression profiles of a large 
number of genes in clinical HCC materials is an essentia) step toward 
clarifying the detailed mechanisms of hepatocarcinogenesis and dis- 
covering target molecules for the development of novel therapeutic 
drugs. 

^ cDNA microarray technology, which enables investigators to ob- 
tain comprehensive data with respect to gene-expression profiles, is 
progressing rapidly. Several studies have already demonstrated the 
usefulness of this technique for identifying novel cancer- related genes 
and for classifying human cancers at the molecular level (4. 5). 
In this paper, we report the identification of genes the expression of 
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which has been ahered during hepatocarcinogenesis through ihe use 
of a genome-wide cDNA microarray containing 23,040 genes. Ex- 
pression profiles of these genes in 20 primary HCCs fell into three 
categories that correlated wel) wiih the infeaion status and type of 
hepatitis virus. AnaJyses of these profiles along with clinicopatholog- 
ical data also facilitated identification of genes associated with tumor 
differentiation and vessel invasiveness. This large body of information 
not only furthers an understanding of the mechanisms of hepatocar- 
cinogenesis but also reveals novel features of known genes and 
identifies additional biological factors involved in liver cancer. 



MATERIALS AND METHODS 

Patients and Tissue Samples. Primary HCCs and corresponding noncan- 
cerous liver tissues were obtained with informed consent from 20 patients who 
underwent hepaieciomy. Patient profiles were obtained from medical records. 
Serologically, 10 cases were hepatitis B surface antigen-positive and 10 cases 
were HCV-positive. No cases wrjih coinfections of HBV and HCV were 
included in this study. Hisiopathological classification was pcrfonned accord- 
ing to the Edmondson grading system: clinical stages were determined accord- 
ing to the Union International Contre Cancer TNM classification. No signifi- 
cant differences were seen between HBV-posiiive and HCV-posttivc status 
with respect to age, sex. grade of differentiation, vessel invasion, or tumor 
stage. 

cDNA Microarrays. We fabricated a "genome-wide" cDNA microauay 
with 23.040 cDNAs selected from the UniGcne database of the National 
Center for Biotechnology Information. The cDNAs were amplified by reverse 
transcription- PCR using poIy(A)-f RNA isolated from various human organs 
as templates; lengths of the amplicons ranged from 200 to 1100 bp without 
rcpciiiivc or poly(A) sequences. The PCR products were spotted in duplicate 
on lype-7 glass slides (Amcrsham) using an Array Spotter Generation III 
(Amersham). Each slide contained 52 housekeeping genes, to normalize the 
signal intensities of the different fluorescent dyes. 

RNA Preparation, Hybridization, and Acquisition of Data. Frozen 
specimens were serially seciioned in lO-^jim slices and stained with H&E to 
define the anaJyzcd regions. To avoid cross-contamination of cancer and 
noncancerous cells, we prepared these two populations by laser- captured 
microdis.section. Total RNA was exiraaed from each population and then 
amplified using Ampliscribe T7 Transcription Kit (Epicentre Technologies). 
The preparation of probes, hybridization, and scanning was performed as 
described previously (6). The fluorescence intensities of Cy5 (noniumor) and 
Cy3 (tumor) for each target spot were adjusted so that ihc mean Cy5 and Cy3 
intensities of 32 housekeeping genes for each slide were equal. 

Validation of Data. To assess the reproducibility of the normalized inten- 
sity ratios, we compared the log2(Cy3:Cy5 intensity ratio) of the 52 house- 
keeping genes between different slide sets. When the difference between 
normalized logarithmic ratios from two experiments was less than I.O, we 
defined the data as reproducible. The reproducibility was more than 90% when 
ihc intensities ofCyS and Cy5 were both above 25.000. 

Classification of 20 HCCs According to Gene Expression Profiles. We 
applied the hierarchical clustering method to both genes and samples. To 
obtain reproducible clusters, we used only selected genes that passed the cutoff 
filler (both Cy3 and Cy5 signals greater than 2b,000 in more than 80% cases 
examined). The analysis was performed using web-available software ("Clus- 
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icr" and *TrccView") wntlcn by M. Eiscn/ Rclore applying the clustering 
algorithm, the fluortscence ratio Tor each spot was first log-transformed; then 
the data for each sample were centered lo remove expen mental biases. 

Identification of Ccncs Responsible for Clinicopalhological Factors. 
We first arranged the relative expression of each gene (Cy3:Cy5 intensity 
ratio) into one of four categories: up-rcgulated (ratio. >2.0). down-regulated 
(ratio. <0.5). unchanged (ratio, between 0.5 and 2.0). and not expressed (or 
slight expression but under the cutoff level for detection). We used these 
categories lo detect changes in expr\:sstuii that were cuiiuiion among samples 
as well as specific to a certain subgroup. To delect differential ly expressed 
genes, we recorded the number of samples in each category within each 
subgroup, for each gene. Then we calculated the U values of Mann-Whitney 
tests, which measured how the sample distributions hciween subgroups over- 
lap. The number or samples within each group is counied and. according to the 
order of the category, the number of overlapped samples ts incorporated into 
ihe U value. A small U shows that the sample disinhution of the two groups 
is clearly separated, e.g.. commonly up-regubtcd in the HBV group and 
down-rcgulaied in the HCV group. We applied a hierarchical clustering algo- 
rithm to all of the selected genes using hamming distance (edit distance). 



* Intcmci addiess: hup://www.nucroarTays.org/sortwarr. 



RESULTS AND DISCUSSION 

Identification of Genes That Were Differenliy Regulated in 

HCCs. To identify genes generally involved in hepatocarcinogcnesis, 
wc compared expression proHles between 20 HCCs and Iheir corre- 
sponding noncancerous liver tissues by means of cDNA microarray. 
We excluded individual data when Cy3 and Cy5 signals were 
<25,000 because data were not reliable for genes giving low signal 
intensities (see *'Maieiia!s and Methods"). When we applied a cutoff 
signalriniensity ratio of cancennoncancer at 2.0 165 genes including 
69 ESTs were selected as t)eing up- regulated in 75% or more of the 
cases examined (Table I). This list of up-regulated genes contained 
MAP4KI as well as MAPK3, suggesting that activation of the MARK 
pathway is a common feature of hepatocarcinogcnesis. Interestingly, 
expression of several genes associated with mitosis, including 
CDC23. rUBGi, CBXl, CKSI, PCTKh PSMD8, CSEJL TTK, and 
PLKI, was commonly increased in cancer cells. As a cell-cycle 
modulator, CDC23 is a known component of the anaphase-promoting 
complex (APC) and leads to meiaphase/anaphase transition through 
I 



CENOMC-WIDE ANALYSIS OF GENE EXPRESSION IN HCCS 
A ,7« Table 2 Commonly down-regulated genes in HCC 
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degradation of cycltn B. TUBGI (-y-tubulin) and CBXI panicipaie in 
centrosome formaiion (7. 8); CKSI and PCTKh encoding cdc2/cdc28 
kinases, arc essential for activation of the anaphasc-promoting com- 
plex. PSMD8 (26S proteasome subunit p3l) is reportedly responsible 
for activation of these kinases (9). Others have reported that CSEJL, 
7TK, and PLKI are associated with formation of the mitotic spindie 
(7, 10) and thai PLKJ can affect the number of centrosomes when 
■exogcnously expressed (II); overcxpression of PLKI has been cor- 
related with poor prognosis in a subset of human cancers (12). Our 
comprehensive expression data for these genes may account for a high 
incidence of chromosomal instability in HCC, and ihey suggest that 
promotion of the mitotic process is generally involved in hepatocar- 
cinogenesis. Therefore, regulation of these miiosis-associaicd genes 
either by chemoiherapcutic agents or by gene delivery might be an 
effective therapeutic strategy for HCCs. 

We also looked for down-rcgulated genes and found 170 (including 
75 ESTs) that were underex pressed in 65% or more of the HCCs 

2 



examined (Table 2) when we applied a cutoff intensity ratio of 
cancer:noncancer at 0.5. The majority of the down-regulaied genes 
encoded hepatocyte-specific gene products {e.g., complement species, 
amyloid, and albumin) and detoxification enzymes (cytochrome 
P-450 and meialloihioncin families), reflecting de-differentiation of 
cancer cells. Regarding retinoid metabolism. LY6E and RBPI, both of 
which appear lo play roles in retinoid-induced differentiation (13, 14) 
were repressed, as was IGFBP3, which also is involved in the retinol- 
mediated inhibition of HCC development (15). Because retinoid is an 
accepted therapy lo encourage differentiation of cells in acute promy- 
clocytic leukemia and is thought to help prevent development of HCC 
(16), reduced expression of these genes may play a crucial role in 
hepaioca rc i nogc ncsi s . 

Wc identified 69 ESTs that were frequently up-regulated and 75 
that were frequently down- regulated, which indicated that a large 
number of genes of unknown function are also involved in hepato- 
carcinogenesis. 
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Pig I. Overan patterns of expression of 4531 
genes across the 20 HCC samples. Red color, ovcr- 
exprrssion m cartcer cells: grem color, underex- 
prcision in cancer cells, black, unchanged expres- 
sion: gray, no expression was dcieCTol (inrcn5ioes 
or borh Cy3 and Cy5 under the coioff value) Crud 
uoted color patterns correspond lo ihc decrees of 
expression changes E^jch ru^. a g^oe. each col- 
ufM. a HCC sample The Jrrulrajfram r/ the 20 
cmrs at the nght of ihc mainji indicates the decier 
of sirnilaniy between tumor samples demonstraiing 
lhai the tumors aie clustered in thrre groups .rrj 
Mur, or fireeny Sample No 123 is a very well 
differentiated tuntor and docs noi appear to belone 
»o any of the clusters The Jendrosrcm at the tZ 
also ,„d,catc$ the degree of sioulariiy among the 
«JI genes examined by expn:ssion patiems 
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Class,ficat.onorHCCsbyG«neF«prevM„„Pronies. We further 
nv«n«a,ed whchcr cl.nicl HCCs c.h^W he clasMf-ed .„.o groups o„ 
he has-s of ,hc,r gene^^x press. on profiles For .his purpose we us^ 

ec'lS'r, . r'"^^ """'^ -'--^ 
we selecled 4.53 1 genes thai passed the cu.off f.I.er (both cy3 and cv5 

^gnals greater than 25.000). Tl^e overall e,pr«SK.n poncn.s across 20 

HCC samples are shown j„ Rg. I. Hk- analyses resulted in the 

tenng of .dcnt.cal genes spotted ..„ different positions into acja 

cent rows, .nd.cat.ng the reliabil.ty of the expression data The 20 

HCCs exajn.,«d fell into three groups, as the !lendrogran, show^ 

o Jwf r'"'' ^•^^P""^*'"'^ f'" ""^ class.ncat,on. wc earned 

ou Spearrnan rank-conelation tests and examined cl.nicopa.holoS^ 

cla s.ficalion. vascular invasion, intrahc-potu: metastasis, and gender 
of .he patients (data not shown). Hc.uever. only the typ^ of hepam^ 

HBV-pos.,.ve and HCV-posiUve HCC ma> result from distinct mech 



Ove Stalus^ To identify genes responsible for the differences betCn 

HBV.pos.t,ve and HCV-positive lumo,.. wc performed 

ney tests and found that 19 known genes amd 21 EVrr^K J 

Significantly different exptession paneril'^^.w^' 

Among the 19 known genes (Fig. 2). seven (CPX2 CYP2EEPHX1 

AKR,C4. FM03. VCTIM. and VCT^BIO) encode key motcuS 

acnvatmg chemotherapeu.ic drugs or detoxifying xeniiotic'tin" 

Most carcinogens are metabolized by Phase I modification enzymes 
ha, generate activated intermediates that are then detoxified by iSS 

EPHxt''Z^F::;xr "'"^ ' ^"^^-"^ '^^'^^s 

.1,- ^ pro-carcinogens to activated 

C?pt^ "T"'^- ----y'"— amine \ activa ^bl^ 

FPHxi.^Slrj'""" hydrocarbons are aaivatcd by 

ZreJsiT f ^"y- --aseJ 

HCV Is ,- H;^r '""^'"^ '""^ ^"^y-"" exclusively in 

HCV-pos..,ve HCCs. which may suggest that their enhanced expres" 



Fig 2. Nineteen known genes of the 40 that 
expressed between HBV-bascd 
and HCV-based HCCs Changes in relative expres- 
sion are presented in graduated color patterns Red 
ovCTcxprcssion; greer,, undcrcxpression: vrllow 
unchanged expression. The number to the left of 
each row is the U value of thr Mann- Whi.ney test 
and the dendrogram indicates the degree of simi- 
lanty between the genes selected. *. the seven 
genes that encode key enzymes for detoxificaion 
of chemotherapeutic drugs or xenobiotic cajcino- 
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Rg 3 Genes the expression of which is rcloicd 10 HCC progression We used colors corresponding 10 relative gene expression as in Rg. 2. Genes related to Kdmondson eradc {A ) 
and 10 vasclar mvos.on (B) Among ,he 321 g.ncs reb.cd ,0 his.olog.cal grade and 151 genes related .0 vascular invasion 128 and 41 noJLd genes are lis^ Blue 
genes that arc associated wiih both va^oiLir invasion and grade of diffeiviitiauon. . F^^.u t y. oiue. 
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Fig 3. Comimied. 
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sion leads to a greater contribution of careinogenic metabolites to the 
mcchan,snr.s of HCV-specific hepatocareim,genesis. 

On the other hand, expression of UCTJA/. UGT2BI0 and CPX2 
was preferentially repressed in HBV-posiiive HCCs WCTIAI was 

leSf !L' °' "^^-P-"'- HCCs e.an,.ned). but expression 
levels of these genes were unchanged in most HCV-positive HCCs In 
accordance with our obsc^ations. Strassburg et al. (21) have shown 
decreased expression of UGTIAI in HCCs as well as in hepatic 
adenomas, .mphcating UGTIAI in an early step of hepatocarcinogen- 

wh!ch f '^'r'*' ^""^^'^ ^^'^^^^ n conjugation reac^on". 
which are frequently related to dctoxiHcation of the active forms of 
carcinogens. GPX2. a major form of glutathione peroxidase in liver 

acSrin Hrr.T*°;!i'""'- g'"ta.h,one peroxidase 

activity in HCCs has been reponed elsewhere (22). Hence reduced 
activities of these enzymes may rcllccl enhanced exposure of hcDa- 
tocytes to activated carcinogens or radicals. Our results suggest that 
decreased expression of detoxification enzymes may be involved 

FunhTr^' "^''^"■^'"^ «f HBV-spcciric hepa.ocarcinogenesis. 

Funhermore. because UGTIAI also catalyzes glucuronidalion of 
SN-38. an active form of irinotecan (23). HBV.posliive HCCs may 
show greater sensitivity to irinotecan than do HCV-positive HCCs 
Different expression patterns among detoxification enzymes should 
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provide information for optimizing the choice and/or .he dosage of 
anticancer drugs for treating IICC patients on a„ individual bJs 

Mr v '^°'"P3""8 expression profiles between HBV-positive 
and HCV-positive HCCs implied that hepatitis viruses affect «pres 
H r^"" "^^ ^ type-specific manner, invoking 

pan ly different mechanisms of carcinogenesis. Consequently, idcnti 
fiction of genes defining virus-type-specific expression profiles 
Zt'eTrns"" " " virus-type-dependen't ..at- 

Identification of Genes Related to HCC Progression. As in the 
multistep model of adenoma-to-carcinoma kequ^enee accept for 
colorectal tumors. HCCs are considen.-d to develop as well-differen- 

s.aSori""'' "'*^"P™87« '° moderately-to-poorly differentiated 
states (24). A companson of expression profiles between well-differ- 
emiated tumors (Edmondson grade I; „ = 7) and modei^tely to poorly 
differentiated tumors (Edtnondson grade II or III; „ = 13- Re 3^) bv 

rpTF.?"^:""'^"""^ "-la'O'alof 321 ge„« (i^ncluding 

1 93 ESTs) that showed different expression patterns between the two 
histologically dtvided groups. In addition .0 the genes encoding l.ver 
specific proteins, they included genes associated with apoptosis and 
T^P^^"!^"%!1!^!^- Apop«osis-reIated genes including TNFSFJO 
TNF5FJ4. CADD34. CFIAR. CLU, CASP6. pHospLidylserine 



GENOME-WIDE ANALYSIS OF GENE EXPRESSION IN HCO 



receptor (25, 26) were preferentially reduced in modcratcly-io- poorly 
differcnuated tumors, implying that a reduced rale of apoptosis is a 
major characteristic of tumor progression. Genes associated with 
immune systems included MAGECh one of the tumor antigens rec- 
ognized by CTLs. whose expression was also repressed only in 
moderaiely-io- poorly differentiated tumors. Reduced expression of 
genes encoding immune targets may confer a growih advantage by 
allowing tumor cells to escape from immune surveillance. 

Furthermore, we investigated expression profiles with respect to 
vascular invasiveness because vascular invasion is a major factor 
aOeciing metastasis and one of the most useful predictive factors of 
prognosis (27). Genes involved in vascular invasion could also rep- 
resent good candidates for new therapeutic targets. We found that 151 
genes (including 1 10 ESTs) were expressed significamly differently 
between noninvasive {n = 8) and invasive (n = 12) tumors (Fig. 35), 
Among the named genes in this category. 19 were associated with 
both vascular invasion and tumor differentiation, but no apoptosis- 
rclated gene was among them; therefore, reduced apoptosis is likely to 
be correlated with tumor de-differeniiaiion and growih, but not with 
va.scular invasion or metastasis. Genes associated with vascular inva- 
sion contamcd ARHC (RhoC)» which was recently reported to play a 
crucial role in metastasis (28). We also found that two other small 
GTra,sc-rcla!ed genes. ARHGAP8 (RhoGAPS) and ARHGEF6, were 
preferentially down-regulated in invasive tumors. Because the small- 
CTPasc Rho family plays important roles in controlling cell motility 
and focal adhesions (29), alterations of their signaling pathways could 
enhance the migratory and invasive capacity of tumor cells and induce 
tumor invasion and metastasis. Although its function is unknown. 
RhoGAPS is thought to inhibit the Rho signaling pathway: hence, 
reduced expression of ARHGAP8 may also result in Rho-mediated 
tumt>r invasion. Our results suggest that controlling the Rho signaling 
pathway cither by reducing expression of ARHC or by inducing 
ARHGAP8 may suppress tumor invasion and subsequent metastasis. 

The genes and their products represented by the numerous ESTs of 
unknown function thai we classified in the same clusters as genes 
associated with apoptosis or immunity may he useful as novel targets 
for drug discovery or rumor markers. Accumulation of data with 
respect to expression profiles of cancer specimens, clinicopathological 
data, sensitivity to treatment, and prognosis will not only help us to 
understand the precise mechanisms of carcinogenesis but also yield 
practical intormation for identifying optimized therapeutic modalities 
and novel therapeutic targets. 
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